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ABSTRACT 

The possible existence of additional long-period planetary-mass objects in the extrasolar planetary 
systems 47 UMa and 14 Her is investigated. We combine all available radial-velocity data on these 
stars, spanning up to 18 years. For the 47 UMa system, we show that while a second planet improves 
the fit to all available data, there is still substantial ambiguity as to the orbital parameters of the 
proposed planetary companion 47 UMa c. We also present new observations which clearly support a 
long-period companion in the 14 Her system. With a period of 6906±70 days, 14 Her c may be in a 
4:1 resonance with the inner planet. We also present revised orbital solutions for 7 previously known 
planets incorporating recent additional data obtained with the 2.7m Harlan J. Smith Telescope at 
McDonald Observatory. 

Subject headings: stars: individual, 47 UMa, 14 Her - stars: planetary systems - extrasolar planets — 
techniques: radial velocities 



1. INTRODUCTION 

The longest-running radial-velocity surve ys are now 
appro achin g tim e baselines of 16-18 years (Butl er et al.l 
rmMlCochra,n et ablll 9971^ . These surveys now achieve 
internal measurement precisions such that the signals 
from long-period giant planets (P^ lOyr) are no\Y enter - 
ing the realm of detectability ijWittenmver et ;ini20nfiD . 
For example, the High Resolution Spectrograph (HRS) 
on the Hobby-Eberly Telesc ope (HET) now has a ve- 
locity precision of 3-4 m s~^ l)Cochran et alTl2004D . the 
Keck HIRES i s achieving 1-2 ni since its 2004 
CCD upgrade iButler et all l20(^. and the HARPS 
instrument has demonstrated precision better than 1 

l l2nnl . 



m s~ IjLovis et al.M200(Tl) . Of particular interest are pu- 
tative multi-planet systems, wherein the residuals of a 
known-planet's orbit show Keplerian periodicity indica- 
tive of a distant giant planet companion in the system. 
Finding systems which contain long-period giant planets 
as well as planets in closer orbits will address important 
questions about the uniqueness of our own Solar System's 
architecture. 

Long-period planet candidates have been reported for 
47 UMa (P = 7.1 yr, Fischer et al. 2002), 55 Cancri 
(P = 12.4 yr, McArthur et al. 2004), HD 217107 (P > 7 
yr, Vogt et al. 2005), and HD 72659 (P = 9.9 yr, Butler 
et al. 2006). Of particular interest are systems which 
contain both short-period and long-period jovian planets. 
Such systems could provide clues to address the question 
of how the processes of planet formation and migration 
can result in both "hot" and "cold" Jupiters in the same 
planetary system. 

In this paper, we present improved fits to the known 
planets in nine systems using all available data, and we 
investigate the possibility of additional long-period ob- 
jects. In §2, we briefly describe the data acquisition and 



fitting procedure. The results are given in §3, including 
our solutions for seven additional planetary systems in 
the McDonald Observatory planet search program. All 
of our radial- velocity measurements for these objects are 
presented in Tables 5-14. In this paper, we use the terms 
"McDonald" to refer to data taken with the McDonald 
Observatory 2.7 m Harlan J. Smith (HJS) Telescope, and 
"HET" to refer to data taken with the 9.2 m Hobby- 
Eberly Telescope. 

2. OBSERVATIONS AND DATA ANALYSIS 

Data obtained from the McDonald Observatory planet 
search program (Endl et al. 2005) are discussed fully in 
iWittenmver et al. 12006.') . Available published data were 
combined with the McDonald data to fi t Keplerian or- 
bits using GaussFit f.Teffer vs et alJlT98^ . GaussFit has 
the ability to allow the offsets between data sets to be 
a free parameter. Parameters of the primary stars con- 
sidered are given in Table 1 ; masses, [Fe/H], and Tg// 
are from lSantos et all l|2flfl4D . and the chromospheric ac- 
tivity index log R'^k derived from Ca II measure- 
ments from the McDonald Observatory spectra. For each 
object, we searched for periodic signals in the residu- 
als to t he known planet's orbit u sing the periodogram 
method llLomblll976l: IScargldll982r . To assess the statis- 
tical significance of those periods, the false alarm prob- 
abilities (FAP) were calculated using the bootstrap ran- 
domization method detailed bv lKiirster et al.l l)1997f) and 
■Endl eTai] ll200l . 
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3. RESULTS AND DISCUSSION 

The Keplerian orbital fits are shown in Figures 1-4, 
and the orbital parameters implied by those fits are given 
in Table 2. In computing the planetary minimum mass 
M sin i and sem i major axis a, the stellar masses derived 
bv lSantos et"an l|2004fl listed in Table 1 were used, with 
adopted uncertainties of 0.05 Mq. 

3.1. ^7 UMa (=HD 95128): Ambiguities concerning the 
second planet 

iButler fc Marcvl l|1996l ^ first reported the 1090-day 
companion to 47 UMa using data from Lick Observa- 
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tory. With additional ve locity measurements over 13 
years, iFischer et alJ ()2002D announced a long-period sec- 
ond planet, 47 UMa c, wit h a period of 2594±90 days 
and a mass of 0.76 Mjup. iNaef et all jlpol) presented 
ELODIE observations of 47 UMa, and noted that the 
second planet was not evident in their data, which were 
fit well with a single Keplerian model. 

We now fit four data sets for 47 UMa: Lick (Fischer et 
al 2002, N=91), ELODIE (Naef et al. 2004, N=44), 2.7m 
HJS telescope (N=35), and the Hobby-Eberly Telescope 
(HEX, N=201). The HEX data, which consist of multi- 
ple exposures per visit, were binned using the weighted 
mean value of the velocities in each visit. The quadra- 
ture sum of the rms about the mean and the mean in- 
ternal error bar was adopted as the error bar of each 
binned point (N=63). The offset between the overlap- 
ping HJS and HET datasets was used to merge them 
into one, which was used in all fits. The one-planet fit 
(Model 1) and the residuals to that fit are shown in Fig- 
ure 1. We emphasize that this fit includes all available 
published data, over a time span of more than 18 years, 
and includes 195 high-precision measurements obtained 
with the HET High-Resolution Spectrograph at 61 inde- 
pendent epochs, which are given in Table 6. The total 
rms about the combined one-planet fit is 10.4 m s~^. 
The rms scatter about the one-planet fit for each of the 
four datasets is: Lick-10.7 m s'^, ELODIE-13.0 m s'^ 
HJS-13.5 m s-\ HET-4.9 m s'^. A periodogram of the 
residuals of all of the data to the 1-planet fit is shown in 
Figure 2. No clear peak rises above the noise level at any 
period between 2 and 10000 days; the total duration of 
the available data is now about 6900 days (18.3 years). 
While a peak is present at about 2 212 days, clos e to th e 
period reported for 47 UMa c bv IFischer et alJ l|2002|) . 
its false-alarm probabihty (FAP) is 27.2%. 

To further explore the possible presence of 47 UMa c, 
we fit all of the datas ets with a twq - planet model fixed 
at the parameters of IFischer et al.l l|2002|) (Model 2), 
and then repeated the fit allowing all parameters to be 
free except for e and uj of the second planet (Model 3), 
which were fixed at 0.005 and 127°, respectively, after 
IFischer et al.l l)2002|) . No models achieved convergence 
with those two parameters free. The fit obtained by 
Model 3 is shown in the right panel of Figure 2. Model 
2 had a reduced chi-square (xt) of 5.81 and an rms of 
12.7 m s^^ about the fit, whereas Model 3 had a xt oi 
2.11 and an rms of 8.4 m s~^. For comparison, the one- 
planet fit (Model 1) had a xt oi 3.23 and an rms of 
10.4 m s^^. Noting that the poor fit of Model 2 was 
largely due to errors in the period of the inner planet, 
we re-did the fits allowing the parameters for the inner 
planet to be free while f ixing those of t he ou ter planet at 
the values reported bv IFischer et all l)2002f) (Model 4). 
The xt of t'^is fit was 2.68, with an rms of 10.2 m s~^. 
These tests are summarized in Table 3. The free two- 
planet (Model 3) fit was the best of the four, in terms 
of both the goodness-of-fit criterion (xt) and the total 
rms scatter about the fits. The rms about the individ- 
ual datasets for this fit is the following: Lick-8.0 m s^^, 
ELODIE-11.1 m s'^, HJS-10.1 m s'^, HET-5.3 m s~^. 
The parameters for the 47 UMa planetary system given 
in Table 2 are those obtained by Model 3. Although the 
best-fit set of parameters obtained a period of 7586 days 
for the outer planet, there is no corresponding peak on 



the periodogram shown in Figure 2. However, the pe- 
riodogram method is not as reliable when the periodic 
signal approaches or exceeds the total duration of obser- 
vations, as is the case here, where the total time baseline 
is 6942 days. We note t hat we have be en able to repro- 
duce the result of Fischer et al.l l)2002t) by this method; 
i.e. a periodogram analysis of the Lick data alone after 
removing 47 UMa b revealed a strong peak at 2083 days, 
but still with a FAP of 0.15% . We are also able to re- 
cover the IFischer et al.l l|2002|) parameters of 47 UMa c 
from the Lick data alone. S ince the total time c overage 
of the Lick data presented in lFischer et alJ l)2002t) is 5114 
days (14 yr), an analysis of those data alone is not ad- 
equate to fully constrain the 7586-day period obtained 
by our best-fit model which adds four years to the total 
duration of observations. It is p ossible that the shor ter 
period for 47 UMa c reported bv IFischer et all lj2002D is 
an alias of the true period; at present, our fits indicate 
that period to be about three times longer, but with sub- 
stantial uncertainty. 

To further test the methods by which we conclude that 
the pa rameters of 47 UMa c reported by IFischer et alJ 
are dubious, we performed some Monte Carlo sim- 
ulations. From each of the three data sets considered in 
the fits described above, we generated 1000 simulated 
sets of velocities consisting of a Keplerian signal plus 
7 m s~^ of Gaussian noise for each of the two planets. 
The parameters for the inner planet were those listed 
in Table 2, and those of the outer planet were those 
of IFischer et alJ l|2002fl . These simulated datasets re- 
tained the times of observation and the error bars of 
the originals. The simulated data were then fit with a 
one-planet model exactly as described above, then the 
residuals of the one-planet fit were examined by the pe- 
riodogram method, to determine whether the signal of 
the second planet was recovered. The criteria for re- 
covery were that the period of the second planet had 
to be detected correctly and with a FAP of less than 
0.1%. Thi s FAP was computed usi ng the analytic FAP 
formula of iHorne fc BaliunasI l|1986D . Of the 1000 trials, 
only 6 did not result in a successful recovery of the signal 
of the second planet. The correct period was recovered 
999 times, and the FAP exceeded 0.1% only 5 times; the 
worst FAP was 0.24%. For comparison, the analytic FAP 
of the 2212-day peak in the residuals of the 1-planet fit 
is 1.7%, a factor of 7 higher. These results indicate that 
our method should have been able to detect the signal 
of 47 U Ma c, had it b e en pr esent with the par ameters 
given bv IFischer et alJ l)2002|) and iButler et alJ £2006). 
We conclude that while an additional long-period object 
may be present, the data currently available do not pro- 
vide sufficient evidence for an orbital solution. 

3.2. 14 Her (~HD 14-5675): Evidence for an outer 
companion 

Thirty-five radial-velocity measurements of 14 Her 
(=HD 145675) obtained at McDonald Observatory 
were combined wit h published d ata from Keck H IRES 
lIRutler et 31.120061) and ELODIE llNa,ef et 31.120041) . The 

fit to the combination of these three data sets is shown 
in Figure 3, and the system parameters implied by that 
fit are given in Table 2. It is evident from Figure 3 that 
the single Keplerian fit is inadequate (rms=13.0 m s^^); 
indeed, the residuals to the fit show a clear curvature. 
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iNaef et aD l)2004|) noted an upward linear velocity 
trend of 3.6±0.3 m s~^ yr~^ in their observations which 
covered the time interval 1994 to 2003, but the Keck data 
of Butler et al. (2003) indicated no such trend. When the 
complete 12 yr duration of observations is examined, we 
see that those Keck data were serendipitously obtained 
during a period of relatively constant residual velocity, 
and the more recent McDonald data are now indicating a 
downward trend. Fitting a double-Keplerian model (Fig- 
ure 4) with To for the outer body fixed at 2449100.0 yields 
an eccentricity consistent with zero: e — 0.02±0.06. Fix- 
ing the eccentricity of 14 Her c at 0.0 gives a minimal 
orbit solution with period P — 6906 ± 70 days, and a 
velocity semi-amplitude K — 24.5 ± 1.4 m s^^. These 
orbital elements imply a minimum mass M sin i=2.1 
Mjup at a semimajor axis of 6.9 AU. This represents 
the lowest possible mass and the shortest period for the 
outer companion. The rms scatter about the three data 
sets is as follows: ELODIE - 10.4 m s"\ Keck - 2.9 
m s~^, McDonald - 5.8 m s~^. The total rms about 
this two-planet fit is 8.4 m s~^, with a xt of 1.67. The 
minimal orbit proposed above for 14 Her c has a perias- 
tron distance of 6.2 AU, and 14 Her b has an apastron 
distance of 3.76 AU. The orbit of 14 Her c, at a large 
semimajor axis with a small eccentricity, is similar to 
that of HD 72659 b ( Butler et al. 2003} and HD 50499b 
l)Vogt et alJl2005( l. We emphasize that such fits are pre- 
liminary, and are only given in order t o plac e a low er limit 
on the period and mass of this object. iCoz dziews ki et aD 
l)200fiD performed dynamical simulations of the 14 Her 
system, and proposed the existe n ce of 14 Her c based 
on fits to data from lButler et al.l 1)20031) and INaef et al.l 
l)2004ft . The two most- favored models had 14 Her c near 
the 3:1 or 6:1 mean-motion resonance (MMR), which co- 
incided with highly stable orbital configurations as de- 
termined by their dynamical simulations. The prelimi- 
nary fit obtained in this work suggests a 4:1 resonance 
for 14 Her c. We attempted to fit a double-Keplerian 
model with 14 Her c fixed at the 3:1 MMR parameters 
indicated in i Gozdziewski et aD l)2006t) . but obtained a 
substantially worse fit {xt = 2.2). 

An upper limit on the mass of 14 Her c can be es- 
timate d from the results of adapt i ve-op tics (AO) imag- 
ing by iLiihman fc Javawardhanal l|2002D , who used the 
Keck II AO system to search for companions to 25 
extrasolar planet-host stars. No candidate compan- 
ions were found around 14 Her, and the detection lim- 
its derived from their study exclude stellar-mass ob- 
jects objects at orbital separations ^ 9 AU (see Luh- 
man & Jayawardhana 2002, their Fig. 10). A simi- 
lar study by [Patience et al.] l)2002() using the Lick AO 
system also excluded stellar-mass objects beyond about 
12.7 AU . We can therefore use t he upper limit pro- 
vided bv lLuhman fc Javawardhanal ((2002) and our lower 
bound to constrain the mass of the outer companion be- 
tween about 2.1 and 80 Mjup. A more definitive state- 
ment on its nature, of course, requires many more years 
of observations, until a second velocity turnaround is con- 
firmed. The large separation implied by these data make 
14 Her an attractive target for future direct imaging at- 
tempts and for astrometric follow-up. The astrometric 
perturbation due to 14 Her c would be a =0.9 mas, us- 
ing the minimal orbit solution given above. For sin i—0.5, 



this sig nal would be a —1.8 mas, equivalent to that of 
e Eri b l)Benedict et al.l l2006'l. Since the astrometric per- 
turbation increases with semimajor axis and planet mass, 
the signal could be much larger (Sozzctti 2005). 

3.3. Revised orbital parameters for 7 additional systems 

The addition of new data from the HJS Telescope pre- 
sented in this paper, and the use of multiple independent 
data sets in fitting Keplerian orbits, have generally im- 
proved the precision of the derived planetary parameters 
by 25-50%. In particular, the precision of the orbital pe- 
riods have been improved by the addition of new data, 
due to the increased number of orbits now observed. In 
this section, we briefly describe the results of our com- 
bined fits for seven additional systems for which the re- 
vised orbital solutions are in agreement with previously 
published results. Our parameters are within 2a of pre- 
vious values except for Tq and u of v And c, and P and 
w of u And d, which differ from lButler e t al. (2006) by 
between 2.7 and 3.8 sigma. 

The inner planet of HD 217107, w ith a period of 7.1 
days, was first reported in lFischer et al. (1999), and ad- 
ditional COR ALIE data supporting this discovery were 
published by INaef et aD (|T001). A shallo w paraboli c 
trend in the residuals was noted bv lFischer etal\ l)2001'l. 
I Vogt et ahl {2005) used more recent Keck data to postu- 
late an outer companion with P > 7 years. However, this 
object has not completed a full o rbit, and hence there is a 
wide range of possible solutions. iWright et al.l l)2005|) re- 
vised the period of HD 217107c to 20,140 days (55 years), 
and Vogt et al. (2005) give a period of 3150 days; our fits 
with the short period improved xt by only 0.07 over the 
long period. We fit four dataset s for HD 217107 : Lick 
( Fischer et a]"rT99 9). CORALIE llNa,ef et al .120011) . Keck 
( Vogt et alJl2005j) . and 20 observations from McDonald. 
We fit a doub le-Keplerian mode l using the parameters for 
planet c from lVogt et alJ l|2005tl as a starting point for the 
least-squares fitting procedure. The rms of this double- 
Keplerian fit is 9.0 m s~^, and the results are given in 
Table 2. The resulting parameters arc in agreement with 
Vogt et al. (2005) to within 2a. Vogt et al. (2005) per- 
formed separate fits to the Lick and Keck data and com- 
mented on the high degree of uncertainty in the parame- 
ters of HD 217107c. The combined fits given i n this work 
support the ~3150-day period reported by I Vogt et all 
(2005), but the uncertainties remain large. At present, 
no definitive statements can be made until this object 
completes a significant fraction of an orbit, or until such 
time as it can be detected via direct imaging or astrome- 
try, techniques which are well-suited to very long-period 
objects. 

The three planets around v And (=HD 9826) were 
fit wit h a triple-Keplerian model, combining the Lick 
data of lButler et all l)2006j) and the Advanced Fiber Op- 
tic Echelle (AFOE) data of lButler etH\ l|1999D with 41 
observations from McDonald. Lick data preceding the 
Hamilton spectro meter upgrade (1995 February) were 
excluded following IB^itler et al.l l|1997D . The total rms 
about the fit is 14.6 m s~^. As shown in Table 4, analy- 
sis of the residuals to a triple-Keplerian fit resulted in 
only a marginally significant periodicity at 2000 days 
( FAP=0.4%). 

iMarcv et all ljl999|) first detected the inner compan- 
ion to HD 168443 using the HIRES spectrograph on 
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Keck I. The outer 4.8-yr planet wa s then re ported in 
iMarcv et all ^200^ and confirmed bv lUdrv et al. ( 2002|. 
A two-planet fit to the M cDo nald data c ombined with 
that of Butler et all l|2006|) and lUdrv et~al . ( 2002) yields 
an rms of 9.5 m s~^. The parame ters given in Table 2 
are in agreement with those of iButler et al.l 1)200(11) within 
2(7. Periodogram analysis of the residuals to the double- 
Keplerian fit revealed a 63-day periodicity with a FAP 
of 0.12% (Table 4). Performing the same analysis on the 
three datasets separately, however, showed no enhanced 
power at that period. Additionally, the period of the in- 
ner planet is 58 days, preventing any object in a 63-day 
orbit. 

The following five planets showed no residual period- 
icities of interest. The hot Jupiter orbiting HD 179949 
was found by the Anglo- Australian Planet Search pro- 
gram using th e 3.9m Anglo-Australian Telescope (AAT) 
l|Tinnev et all 12001(1 . The rms about the combined 
fit is 11.5 m s~^, and the fitte d parameters agree 
with those of IB hitler et al.l l)2006 ^ to within la. For 
16 Cyg B (=HD 1864 27), data from the dis covery paper 
llCochran et alJ Il997j) and Lick data from iButler et al.l 
l(2006tl were combined with 37 additional measurements 
from McDonald Observatory. The rms about the com- 
bined fit is 10.6 m s~^, and the orbital elements are 
within 2a of those given in Butler et al. ( 2006) . For 
HD 195019b, Lick data from iButler et al] ' l(200(l were 
combined with 19 McDonald observations. The rms 
about the combined fit is 15.8 m s~^ , and the param- 
eters agree with those of (Butler ct al.' f2006') within 2a. 
For HD 210277b, data from Butler et al. (2006) and 
INaef et al.l Emi) were combined with 21 measurements 
from McDonald. The rms about the combined fit is 6.8 
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4. SUMMARY 

We have combined new radial-velocity observations 
from McDonald Observatory with previously published 
data to improve the precision of the orbital parameters 
for 9 planetary systems. The largest data set yet as- 
sembled on 47 UMa indicates no statistically significant 
signal attributable to a second planet in the system. We 
have examined the residuals to our Keplerian fits, and for 
the case of 14 Her, we find clear evidence for a distant 
outer companion, which appears to be in a 4:1 resonance 
with the inner planet. 
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TABLE 1 



Stellar Parameters 


Star 


Spec. Type 


Mass (Mq) 


[Fc/H] 


Teff (K) 




V And 


F8V 


1.30 


0.13±0.08 


6212±64 


-5.01 


47 UMa 


GIV 


1.07 


0.06±0.03 


5954±25 


-5.04 


14 Her 


KOV 


0.90 


0.43±0.08 


5311±87 


-5.06 


HD 168443 


G5 


0.96 


0.06±0.05 


5617±35 


-5.14 


HD 179949 


F8V 


1.28 


0.22±0.05 


6260±43 


-4.75 


16 Cyg B 


G3V 


0.99 


0.08±0.04 


5772±25 


-5.03 


HD 195019 


G3IV-V 


1.06 


0.08±0.04 


5842±13 


-4.89 


HD 210277 


GOV 


0.92 


0.19±0.04 


5532±14 


-5.11 


HD 217107 


G8IV 


1.02 


0.37±0.05 


5646±26 


-5.13 



TABLE 2 
Keplerian Orbital Solutions 



Planet 


Period 


To 




e 




K 




M sin i 




a 






(days) 


(JD-2400000) 






(degrees) 


(m 


s~i) 




(Mj,p) 




(AU) 




V And b 


4.61708±0.00006 


50001. 8±0.4 





029±0.013 


46±29 


71 


1±1 





0.69±0 


03 


0.059±0 


001 


V And c 


241.52±0.21 


50149. 7±3.3 





254±0.016 


232.4±4.9 


56 


1±1 


2 


1.98±0 


09 


0.83±0 


01 


V And d 


1274.6±5.0 


50074±16 





242±0.017 


258.5±5.4 


64 


1±1 


1 


3.95±0 


16 


2.51±0 


04 


47 UMa b 


1083.2±1.8 


50173±65 





049±0.014 


111±22 


49 


3±1 





2.60±0 


13 


2.11±0 


04 


47 UMa c'' 


7586±727 


52134±146 





005 (fixed) 


127 (fixed) 


13 


3±1 


4 


1.34±0 


22 


7.73±0 


58 


14 Her b 


1773.4±2.5 


51372. 7±3.6 





369±0.005 


22.6±0.9 


90 


0±0 


5 


4.64±0 


19 


2.77±0 


05 


HD 168443 b 


58.1112±0.0009 


50047.45±0.04 





530±0.001 


172.7±0.2 


476 


Oil 





7.48±0 


27 


0.290±0 


005 


HD 168443 c 


1765.8±2.2 


50255±4 





222±0.003 


64.6±0.8 


299 


2±1 





16.87±0 


64 


2.84±0 


05 


HD 179949 b 


3.09250±0.00003 


51793.98±0.29 





022±0.014 


183±34 


112 


8±1 


6 


0.95±0 


04 


0.045±0 


001 


16 Cyg B b 


799.5±0.6 


50539. 3±1.6 





689±0.011 


83.4±2.1 


51 


2±1 


1 


1.68±0 


07 


1.68±0 


03 


HD 195019 b 


18.2008±0.0003 


50033. 1±0.8 





014±0.004 


239±16 


272 


8±1 





3.67±0 


13 


0.138±0 


002 


HD 210277 b 


442.1±0.4 


50988. 2±1.6 





472±0.011 


118.2±1.9 


39 


5±0 


5 


1.23±0 


05 


l.lOiO 


02 


HD 217107 b 


7.12689±0.00005 


49998.50±0.04 





132±0.005 


22.7±2.0 


140 


6±0 


7 


1.33±0 


05 


0.073±0 


001 


HD 217107 c 


3352±157 


50921±84 





537±0.026 


164 (fixed) 


39 


8±6 


4 


2.50±0 


48 


4.41±0 


21 



''Results for 47 UMa from Model 3 in Table 3. 



TABLE 3 
47 UMa Orbital Solutions 



Parameter 


Model 1 


Model 2=^ 


Model 3*^ 


Model 4= 


Pb (days) 


1078.3 


1089.0 


1083.2 


1073.1 


Tj, (JD-2400000) 


52391 


50356 


50173 


50148 




0.088 


0.061 


0.049 


0.028 


u>i, (degrees) 


127 


171.8 


111 


100 


Kt (m s-i) 


46.7 


49.3 


49.3 


50.2 


Pc (days) 




2594 


7586 


2594 


Tc (JD-2400000) 




51363.5 


52134 


51363.5 


ec 




0.005 


0.005 


0.005 


LJc (degrees) 




127 


127 


127 


(m s-i) 




11.1 


13.3 


11.1 


rms (m s^-'-) 


10.4 


12.7 


8.4 


10.2 




3.23 


5.81 


2.11 


2.68 



^All parameters fixed at those of lFischer et all 120021) 

'^All parameters free except for Cc and uic ^^^^^^^^^^ ^^^^ 

'^Oiily parameters for planet c fixed at those of lFischer et all i2002H 
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TABLE 4 
Periodogram Analysis 



Star 


Period (days) 


FApa 


V And 


2000.0 


0.004 


47 UMa 


2212 


0.272 


14 Her'' 


2.37 


0.118 


HD 168443 


63.37 


0.001 


HD 179949 


32.05 


0.023 


16 Cyg B 


2.32 


0.895 


HD 195019 


41.60 


0.045 


HD 210277 


15.50 


0.392 


HD 217107'' 


485.44 


0.332 



''10000 bootstraps. 

''Residuals obtained from 2-planet fit. 



TABLE 5 

HJS Telescope Radial Velocities for v And 

(-IID 9826) 



JD-2400000 Velocity (m s'^) Uncertainty (m s"') 



51452.86152 


-13.8 


9.2 


51452.86540 


-3.0 


13.4 


51504.69248 


0.8 


8.1 


51530.78235 


-22.2 


9.4 


51532.65475 


-16.0 


10.8 


51557.64589 


26.8 


7.7 


51750.96883 


113.1 


10.6 


51775.91401 


-54.7 


8.1 


51778.93251 


63.8 


9.8 


51809.77554 


-62.0 


8.6 


51859.76793 


-50.3 


6.8 


51861.85297 


96.7 


7.1 


51861.85893 


96.1 


7.7 


51862.79111 


91.5 


11.9 


51862.79735 


68.8 


10.1 


51918.75978 


5.0 


8.7 


51918.76434 


23.1 


9.5 


51920.67309 


46.5 


6.9 


51946.68903 


23.8 


14.3 


51987.57597 


6.3 


9.1 


52142.87690 


55.5 


9.3 


52218.75362 


-63.5 


9.6 


52249.68686 


-6.3 


8.5 


52329.57235 


-125.9 


8.7 


52495.93509 


-151.8 


14.9 


52539.84157 


-104.7 


9.3 


52539.84425 


-113.3 


9.4 


52577.85533 


5.0 


8.5 


52619.75862 


44.4 


8.5 


52933.79237 


102.3 


10.4 


52933.79503 


92.7 


9.8 


52958.68304 


-5.6 


9.1 


53015.67064 


28.9 


8.2 


53035.57199 


-20.0 


7.0 


53394.63389 


77.3 


10.2 


53394.64950 


85.3 


10.0 


53394.65545 


76.5 


10.8 


53632.91257 


-106.3 


9.2 


53690.78040 


-29.4 


10.0 


53691.77259 


-130.3 


13.4 


53746.71222 


-150.9 


8.3 
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TABLE 6 

HJS Telescope Radial Velocities for 47 UMa 

f-IID 95128) 



JD-2400000 Velocity (m s"!) Uncertainty (m s'^) 



51010.62898 


51.6 


6.4 


51212.97474 


-10.8 


5.5 


51240.81250 


-7.3 


6.2 


51274.78993 


-11.2 


5.4 


51326.70558 


-22.4 


6.2 


51504.95996 


-42.2 


6.0 


51530.01978 


-29.4 


6.9 


51555.94972 


-25.5 


5.8 


51655.74023 


5.8 


5.8 


51686.75156 


-6.7 


6.5 


51750.60418 


2.0 


6.6 


51861.01895 


53.7 


6.9 


51917.93086 


47.5 


7.1 


51987.85527 


47.8 


8.5 


52004.83235 


59.7 


6.0 


52039.77936 


54.8 


7.5 


52116.60554 


39.4 


7.6 


52249.00010 


9.5 


7.6 


52303.89238 


-9.6 


5.5 


52305.84757 


-11.8 


6.1 


52327.86285 


12.3 


16.6 


52353.85949 


-12.9 


7.7 


52661.95399 


-24.9 


5.4 


53017.93695 


61.5 


7.5 


53069.76686 


60.7 


6.4 


53692.03243 


-49.7 


8.1 


53748.89147 


-47.9 


6.0 


53787.91198 


-35.2 


6.3 


53805.88756 


-29.3 


5.6 


53809.80777 


-30.8 


6.1 


53805.88756 


-29.3 


5.6 


53809.80777 


-30.8 


6.1 


53787.91198 


-35.2 


6.3 


53861.74397 


-17.1 


6.1 


53909.61977 


13.8 


7.0 
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TABLE 7 

HET Radial Velocities for 47 UMa (=HD 95128) 



JD-2400000 


Velocity (m s ^ ) 


Uncertainty (m s ^) 


53313.99225 


61.3 


1.5 


53313.99417 


64.1 


1.6 


53313.99608 


52.5 


1.7 


53314.99012 


56.2 


1.3 


53314.99203 


59.9 


1.3 


53317.98891 


48.4 


1.6 


53317.99082 


48.4 


1.4 


53317.99273 


47.5 


1.4 


53334.94873 


48.2 


1.4 


53334.95296 


50.9 


1.4 


53334.95591 


54.0 


1.4 


53335.94413 


52.8 


1.4 


53335.94708 


54.5 


1.4 


53335.95003 


54.5 


1.3 


53338.92569 


57.9 


2.0 


53338.92756 


53.0 


3.1 


53338.92947 


46.4 


2.1 


53338.93800 


45.0 


1.9 


53338.93991 


48.6 


1.9 


53338.94181 


57.3 


1.8 


53338.94426 


52.2 


1.9 


53338.94617 


46.5 


1.8 


53338.94808 


47.7 


1.9 


53340.91533 


58.2 


1.4 


53340.91724 


53.0 


1.5 


53340.91915 


58.0 


1.4 


53346.92015 


50.0 


1.5 


53346.92207 


50.6 


1.4 


53346.92398 


51.7 


1.3 


53348.90749 


51.3 


1.6 


53348.90939 


56.8 


1.6 


53348.91131 


52.4 


1.5 


53350.91699 


51.9 


1.4 


53357.87818 


39.1 


1.8 


53357.88009 


50.8 


1.8 


53357.88200 


40.8 


1.6 


53359.87351 


46.6 


2.2 


53359.87542 


43.8 


2.1 


53359.87733 


48.4 


2.1 


53365.86302 


39.9 


1.9 


53365.86489 


48.7 


2.1 


53365.86678 


47.2 


2.0 


53367.86198 


49.4 


1.8 


53367.86389 


41.8 


1.7 


53367.86580 


43.7 


1.8 


53371.85542 


32.8 


1.8 


53371.85733 


41.6 


2.5 


53371.85925 


41.8 


1.7 


53373.85759 


40.9 


2.1 


53373.85950 


32.0 


4.0 


53389.79570 


42.8 


2.0 


53389.79762 


37.9 


2.0 


53389.79953 


35.2 


1.9 


53391.79094 


35.3 


2.0 


53391.79285 


36.6 


2.1 


53391.79477 


34.5 


1.8 


53395.77629 


42.2 


1.8 


53395.77819 


32.7 


1.8 


53395.78010 


42.7 


1.9 


53400.99279 


33.2 


1.7 


53400.99470 


28.6 


1.7 


53400.99661 


26.1 


1.7 


53408.76776 


37.1 


1.9 


53408.76968 


30.5 


2.1 


53408.77158 


29.1 


2.2 


53414.72643 


28.5 


2.1 


53414.72832 


34.6 


2.0 


53414.73023 


33.9 


1.9 


53416.70849 


29.2 


2.2 


53416.71038 


31.4 


2.1 


53416.71231 


33.6 


2.1 


53421.93924 


23.3 


1.3 


53421.94115 


25.9 


1.4 


53421.94306 


24.6 


1.4 


53423.70290 


29.4 


1.3 


53423.70481 


26.5 


1.4 


53423.70672 


23.0 


1.4 


53432.90612 


21.4 


1.3 


53432.90802 


18.9 


1.3 
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TABLE 8 

HJS Telescope Radial Velocities for 14 Her 

(-IID 11507.-,) 



JD-2400000 


Velocity (m s ^) 


Uncertainty (m s ^) 


51329.82953 


154.3 


7.4 


51358.76587 


148.8 


6.5 


51417.72032 


127.8 


7.5 


51449. o2d44 


106.7 


6.9 


51656.89886 


11.7 


7.3 


51689.69367 


-0.5 


7.2 


51750.74058 


-0.8 


8.0 


51777.71797 


-20.2 


6.7 


51809.60771 


-21.1 


6.8 


52037.88607 


-31.4 


7.0 


52115.73144 


-35.1 


7.1 


52145.70159 


-23.4 


7.7 


ozioi.Dzyio 


-23.7 


7.1 


52354.98989 


-25.0 


12.7 


52451.80971 


-2.1 


7.4 


52454.80637 


-3.2 


6.7 


52471.72725 


-12.3 


6.4 


52495.71816 


-13.9 


7.4 


52541.58879 


-0.7 


7.0 


52806.74254 


53.6 


6.9 


52841.82135 


55.6 


7.5 


52933.55588 


90.8 


7.2 


53215.73683 


90.9 


8.0 


53505.84769 


-31.7 


10.0 


53505.87913 


-42.9 


8.2 


53563.74131 


-39.8 


6.6 


53586.66183 


-31.3 


8.0 


53633.61295 


-47.0 


7.2 


53636.60040 


-56.8 


6.6 


53805.95471 


-71.0 


7.2 


53809.96897 


-66.3 


7.2 


53840.90354 


-62.0 


8.2 


53863.84247 


-62.0 


9.3 


53909.73446 


-57.6 


7.3 


53927.75413 


-58.5 


6.8 



TABLE 9 

HJS Telescope Radial Velocities for HD 168443 



JD-2400000 Velocity (m s"i) Uncertainty (m s^^) 



51329.88720 


-523.3 


6.2 


51360.83723 


161.3 


5.9 


51417.74908 


208.1 


6.2 


51451.63634 


-54.8 


6.2 


51689.83601 


267.1 


5.9 


51752.70677 


369.2 


5.8 


51776.68347 


364.6 


6.7 


51810.65202 


369.8 


5.8 


51861.53311 


252.0 


7.5 


52040.91217 


136.4 


5.7 


52116.78446 


97.1 


6.5 


52453.83556 


-197.4 


6.1 


52473.72399 


-233.8 


6.4 


52492.69696 


-766.5 


7.6 


52540.69172 


-552.2 


6.3 


52577.57953 


-156.6 


5.8 


52840.83905 


-684.6 


7.1 


52932.59046 


2.3 


5.7 


53504.89012 


392.3 


6.5 


53565.82109 


442.8 


6.8 


53863.92552 


78.8 


6.8 


53911.89301 


27.4 


6.6 
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TABLE 10 

HJS Telescope Radial Velocities for HD 179949 



JD-2400000 


Velocity (m s~ ) 


Uncertainty (m s~ ) 


O J-UU 1 .\) 1 uto 


-72.4 


7.7 


■S1 121 ^'^727 


109.3 


10.2 




94.1 


8.3 


"^1 '^fiO SI 41 Q 


-106.4 


9.0 




oo.u 


^ 9 
o.z 




72.1 


9.0 


51750.78249 


-121.0 


9.2 


51775.71719 


-81.5 


8.6 


51812.66687 


-89.6 


11.2 


52040.93486 


-62.9 


8.6 


52116.80005 


104.0 


9.2 


52492.71373 


-77.8 


8.7 


52577.51315 


98.7 


8.4 


52840.82950 


100.4 


12.2 


52933.58832 


122.3 


8.0 


53565.85370 


-104.1 


10.6 


53911.90672 


-41.0 


10.4 



TABLE 11 

HJS Telescope Radial Velocities for 16 Cyg B 

f-IID 1R6127) 



JD-2400000 


Velocity (m s ^) 


Uncertainty (m s ^) 


51008.86449 


-7.0 


6.7 


51011.86885 


-1.9 


7.0 


51065.73704 


0.5 


7.7 


51121.69528 


15.6 


7.9 


51154.57100 


14.7 


7.3 


51328.94395 


29.1 


6.7 


51361.87182 


-43.2 


7.5 


51414.67124 


-48.4 


7.3 


51449.78342 


-41.3 


6.8 


51502.60158 


-42.5 


7.8 


51529.57703 


-36.6 


8.0 


51689.90882 


-19.1 


6.9 


51753.74228 


7.5 


7.6 


51777.79311 


-3.6 


6.9 


51861.66676 


24.4 


7.1 


520.39.91285 


35.7 


6.4 


52115.84910 


48.7 


6.9 


52144.77541 


-16.5 


7.0 


52181.72467 


-51.1 


7.3 


52451.88739 


-16.3 


6.2 


52471.85446 


-15.1 


6.4 


52495.78667 


2.4 


9.2 


52538.66450 


-7.3 


6.8 


52577.67137 


-9.1 


7.4 


52600.58553 


1.5 


7.4 


52620.56341 


0.2 


7.4 


52807.88580 


23.2 


6.8 


52840.96694 


34.9 


7.7 


52930.68220 


36.5 


7.5 


52932.66273 


21.3 


6.7 


52960.59294 


-27.4 


7.3 


53163.86157 


-7.0 


7.0 


53321.64999 


-9.6 


14.7 


53585.84777 


27.6 


7.0 


53654.68515 


45.4 


7.0 


53689.59980 


57.8 


7.4 


53907.87292 


-24.1 


7.7 
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TABLE 12 

HJS Telescope Radial Velocities for HD 195019 



JD-2400000 


Velocity (m s~ ) 


Uncertainty (m s~ ) 




-165.6 


5.6 




-199.1 


6.2 


O i- i 1 \} . i \}0:J0 


-221.6 


4.8 


77^ 747'^7 


-191.2 


4.5 


O-LoUU-UlOZU 






O-LOUZr.UliJOO 


152.9 


6.3 


^21 1 4 Qnn'S2 

rj^ ±. J_rt» V J O ^ 


317.7 


5.9 




290.2 


7.1 


52472.85791 


-14.4 


6.9 


52492.82837 


136.8 


7.5 


52538.72733 


-74.6 


6.4 


52807.86835 


268.4 


6.9 


53321.56421 


-95.5 


12.3 


53505.93331 


-185.1 


22.4 


53563.87408 


-103.1 


5.6 


53584.81430 


159.5 


7.8 


53633.65906 


-217.5 


7.7 


53636.72312 


-86.5 


6.7 


53691.64172 


-69.9 


6.8 



TABLE 13 

HJS Telescope Radial Velocities for HD 210277 



JD-2400000 Velocity (m s'^) Uncertainty (m s"!) 



51531.59673 


-30.4 


5.4 


51557.53626 


13.1 


6.1 


51558.54314 


-1.3 


6.3 


51689.93883 


12.3 


6.3 


51750.88078 


24.3 


5.9 


51776.81672 


27.7 


5.3 


51860.63668 


1.9 


11.7 


51917.54532 


-58.3 


7.4 


52116.88917 


10.3 


7.7 


52221.67895 


30.4 


6.2 


52473.83462 


-9.8 


5.6 


52492.85694 


-2.7 


8.2 


52540.77428 


6.2 


6.4 


52621.61830 


9.1 


5.3 


52840.95838 


-28.5 


6.0 


52931.71553 


3.2 


6.1 


53563.94438 


35.3 


5.5 


53633.72641 


-26.7 


5.7 


53635.81499 


-11.6 


5.9 


53689.64535 


-46.1 


5.8 


53927.90889 


41.6 


8.1 



Long-Period Companions 



13 



TABLE 14 

HJS Telescope Radial Velocities for HD 217107 



JD-2400000 


Velocity (m s~ ) 


Uncertainty (m s~ ) 




-44.9 


4.8 


ij i-ij'jyj.'j 1 iyijyj 


-52.7 


5.1 




154 4 


^ 9 


"^1 77f\ S9S(SR 


-79.6 


4 9 


^^ 777 s'^7^'^ 

o± III .oo I yjo 


0.8 


5 7 


^^ 77R ^CiQ^d 


97.5 


4.6 


^1 77S 8227'^ 

iJ-L 1 1 O. Oiiii 1 iJ 


101.5 


5.3 


o±ouy.ouoo± 


38.2 


5.0 




-28.0 


16.0 




-75.8 


6.8 




-102.3 


4.7 


_L ± 1 7 . O t7 .y O O 


-15.7 


7.3 


52219.68045 


1.5 


6.5 


52473.84637 


-83.3 


5.0 


52492.88220 


157.8 


8.7 


52540.78985 


46.9 


5.3 


52895.85865 


-72.5 


5.0 


52931.80527 


-49.1 


4.9 


53017.56950 


-29.9 


6.0 


53563.89812 


-51.8 


4.8 


53630.81136 


-9.2 


6.8 


53635.82600 


-81.0 


6.5 


53689.65961 


177.3 


5.9 



14 



Wittenmyer, Endl, Cochran 





Fig. 2. — Left panel: Lomb-Scargle periodogram for 47 UMa, after removal of 47 UMa b. Right panel: Best-fit double Keplerian model 
for 47 UMa. The symbols have the same meaning as in Fig. 1. 
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Fig. 3. — Same as Fig. 1, but for a one-planet fit for 14 Her b. Open circles are from iNaef et all 120041^ . open triangles are from 
IButler et all [200^ . and filled circles are from McDonald. 
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Fig. 4. — Left panel: Double-Keplerian fi t for 14 Her. Right panel: Residuals of the 2-planet fit. Open circles are from lNaef eraLl Eoti^ , 
open triangles are from IButler et all I200fii) . and filled circles are from McDonald. 



